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There has been much debate over whether or not one could explain the observed acceleration of 
the Universe with inhomogeneous cosmological models, such as the spherically-symmetric Lemaitre- 
Tolman-Bondi (LTB) models. It has been claimed that the central observer in these models can 
observe a local acceleration, which would contradict general theorems. We resolve the contradiction 
by noting that many of the models that have been explored contain a weak singularity at the location 
of the observer which makes them unphysical. In the absence of this singularity, we show that 
LTB models must have a positive central deceleration parameter go, in agreement with the general 
theorems. We also show that it is possible to achieve a negative apparent deceleration parameter 
at nonzero redshifts in LTB models that do not contain this singularity. However, we find other 
singularities that tend to arise in LTB models when attempting to match luminosity distance data, 
and these generally limit the range of redshifts for which these models can mimic observations of 
an accelerating Universe. Exceptional models do exist that can extend to arbitrarily large redshift 
without encountering these pathologies, and we show how these may be constructed. These special 
models exhibit regions with negative effective equation of state parameter, which may fall below 
negative one, but we have failed to find any singularity-free models that agree with observations. 
Moreover, models based on dust-filled LTB metrics probably fail to reproduce observed properties 
of large scale structure. 

PACS numbers: 98.80.-k, 98.80.Jk, 98.80.Es 



I. INTRODUCTION 



The Universe appears to be expanding at an accelerating rate. This has been deduced from luminosity distance 
measurements of Type la supernovae, which appear dimmer than one would expect based on general relativity 
without a cosmological constant 0,0, and from measurements of the current matter density Sl M « 0.27, which is 
too small to close the Universe as required by cosmic microwave background radiation (CMBR) observations 0. 
Many explanations for this discrepancy have been put forward, most entailing a modification of general relativity on 
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cosmological scales or the addition of a new field with exotic properties, called "dark energy"; for reviews, see 

However, there have also been attempts to explain this seemingly anomalous cosmic acceleration as a consequence 
of inhomogeneity rather than modified gravity or dark energy. Many recent papers have claimed that superhorizon 
density perturbations, left over from inflation, could backreact and drive the acceleration of our local universe |a,l3- 
Problems with this argument have been pointed out |8U9tll0 Ull| . and it does not appear to be a plausible alternative. 

In this paper we will only consider subhorizon perturbations, which have also been considered as a possible way 
to obtain an accelerating universe. Kolb et al. use second order perturbation theory 0, 0] to suggest that density 
perturbations could backreact to cause accelerated expansion without the need to introduce any form of dark energy, 
which is an appealing prospect. Rasanen 0] and Notari 0] had earlier explored this claim by looking at average 
expansion parameters in order to argue that we could measure acceleration due to this effect. In opposition to 
this, Siegel and Fry 0] and Ishibashi and Wald 01 have argued that our universe is very accurately described 
with a Newtonianly perturbed Friedmann-Robertson- Walker (FRW) metric, and, treated as such, it does not permit 
accelerated expansion due to perturbations. 

Since it has proven to be quite complicated to analyze the full three-dimensional backreaction problem analyzed by 
Kolb et al., a useful class of models to explore are the spherically-symmetric, yet inhomogeneous, Lemaitre-Tolman- 
Bondi (LTB) 18] cosmological models, containing only cold dark matter, or "dust", and wherein it is often, but not 
always, assumed that we live at the symmetry center. In this way, we can confront the simpler and more general 
question: Are there any models based on general relativity and cold dark matter which can match the observations? 
We cannot completely adress this question with LTB models, which are unrealistic since they place us near the center 
of the Universe, but these models are nevertheless useful toy models to address this general question. More specifically, 
in the LTB models we can ask if a centrally located observer can mistakenly interpret astronomical observations of 
redshifts and luminosity distances as requiring acceleration of the expansion of the Universe. We find that the answer 
is "yes" , and this implies that the mechanism studied by Kolb et al. is somewhat more plausible and requires more 
study. 
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Other papers have used LTB models in analyzing whether or not subhorizon perturbations could backreact and 
drive accelerated exp ansion. Nambu et al. take averages to find effective expansion parameters of specific illustrative 
example models |l9j . Moffat looks at examples [2*3, Mansouri constructs a model that consists of a local LTB patch 
which is embedded into a background FRW spacetime [2ll , and Chuang et al. numerically produce examples of LTB 
models with apparent acceleration [22(. Alnes et al. |2Jj argue against acceleration, but only by looking at a class of 
example models. 

It has been claimed that it is possible to find LTB cosmological models that have qo < 0, where qo = q(z = 0) 
is the deceleration parameter measured by the central observer. However, there are general theorems that prohibit 
such behavior 0, [l3- In Section II, we will first give a general review of LTB models and then we will discuss 
this contradiction and its resolution: there is a local singularity at the symmetry center of models with q < 0, 
corresponding to a non-vanishing radial central density gradient and divergent second derivatives of the density. We 
will prove that excluding this singularity will necessarily lead to a positive value for qo. This singularity is not 
taken into account in any of the above papers, and most of them look at models which are singular at the center 
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We will also show that it is possible to construct models without a central singularity in which one would measure 
negative deceleration parameters q(z), and therefore would measure regions of acceleration, at nonzero redshifts z. 
We will do this by choosing the LTB model and computing the resulting luminosity distance and ultimately q(z); 
we call this the "forward problem". As we discuss in Section II below, LTB models are characterized by two free 
functions of radius r, a bang time function to{r) and an energy function E(r). We focus on LTB models with zero 
energy functions but non-zero bang time functions, because we do not expect the former to produce acceleration. 
This is because, as we will show, the energy function is associated with the growing mode of linear theory, whereas 
the bang time function is associated with the shrinking mode. 

In Section III, we will explore the "inverse problem" , where one chooses the luminosity distance as a function of 
redshift and then attempts to find a corresponding LTB model, which may or may not exist. Here, too, we only 
consider E(r) = 0. We show that there are numerous pitfalls to this method, as other singular behaviors arise which 
generally limit the range of redshifts for which this class of models could reproduce the observed supernova data. For a 
given luminosity distance Dl(z) = rpRw(z)(l + z), there is a critical redshift z crrt where dlnrpRwiz) / 'dln(l + z) = 1. 
For almost all choices of Dl(z), any attempt to find a corresponding zero energy LTB model will fail at some redshift 
smaller than z cr ;t when a singularity is encountered. There are exceptions which pass through a "critical point" at 
z = z cr it, the simple FRW model being one obvious example of such a "transcritical" solution. We show how others 
may be constructed. These models show redshift domains with enhanced deceleration as well as acceleration, but do 
not appear to be consistent with observational data on Dl{z). 

Several papers have already computed how the dependence of the luminosity distance on redshift is distorted in LTB 
models due to purely radial inhomogeneities, and have claimed that we could be tricked into thinking that we are in 
a homogenous accelerating universe when we are really in a dust-dominated inhomogeneous universe |24l l25l l2(i |27| . 
However, this claim has not until now been correctly justified, since all previous papers neglected the central singularity 
and the critical point. 



II. THE FORWARD PROBLEM 
A. Lemaitre-Tolman-Bondi Models 

Using the notation of Celerier [2(|, the LTB spacetime 01 has the line element 

ds 2 = -dt 2 + J^2$-dr 2 + R 2 (r, t) (d8 2 + sin 2 Odcf 2 ) (2.1) 
1 + 2E(r) 

where primes denote derivatives with respect to the radial coordinate r, and E(r) is a free function, called the "energy 
function". We define the function k(r) by k(r) = —2E(r)/r 2 . If k(r) = 0, the Einstein equations admit the solution 

R(r,t) = (67rGp) 1/3 r[t~to(r)} 2/3 , (2.2) 

where fo(r) is another free function, often referred to as the "bang time" function, and p is a fixed parameter. If 
k(r) < for all r, we have the parametric solution 

4nGpr 
R = ] 

t-t (r) = 3[ _ 4 ^, /2 (sinh.-,) , (2.3) 
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and if k(r) > for all r, we have the solution 



AirGpr 
R = — (1 — cosu) 



*-*o(r) 



3fc(r) 
47rGp 

3[*(r)] 



3/2 



(it — sin it) 



(2.4) 



These are Eqs. (18), (19), and (20) of Celerier [2(|, but specialized to the choice M(r) = 47rr 3 /5/3 by choosing the 
radius coordinate appropriately, where M(r) is the mass function used by Celerier, and where p is a constant. 1 The 
energy density of the matter in these models is given by 



P(r,t) 



~ 2 

pr 



R'R? ' 

We define po(t) = p(0,i) to be the central density, and from Eqs. I|2.2(l - (|2.4(l we find 

1 



Po(t) = 



6ttG [t-t o (0)Y 



(2.5) 



(2.6) 



Throughout this paper we will restrict attention to an observer located at r = and at t = t a , where t a is the 
observation time, not to be confused with the bang time to(r). We also choose units such that p = po(t ), and we 
choose the origin of time such that t (0) = 0. A light ray directed radially inward follows the null geodesic 



dt 



R'(r,t) 
\J\ — k(r)r 



-.dr 



and has a redshift given by 



dz = (1 , R'[r,t(r)} 
dr y/l - k(r)r 2 



(2.7) 



(2.8) 



where overdots denote partial derivatives with respect to time and where t[r) is evaluated along light rays that are 
moving radially inward according to Eq. I|2.7|) . Equations 12. 5|) and l|2.8|l give us two important restrictions on the 
derivatives of R(r,t): (i) in order for the density to remain finite, we require R' > 0, which excludes shell-crossing, 
and (ii) in order to have a monotonically increasing z(r), we require R' > 0. 

The luminosity distance measured by the observer at r — and at t = t is given by [2(j 



D L {z) = (l + zf R, 



(2.9) 



where z and R are evaluated along the radially-inward moving light r ay. It is not obvious how to define the deceleration 
parameter in an inhomogeneous cosmology, and Hirata and Seljak [10] explore several definitions. In this paper, we 
restrict our attention to the deceleration parameter that would be obtained from measurements of luminosity distances 
and redshifts assuming a spatially flat FRW cosmology. 2 We can deduce the effective Hubble expansion rate H(z) 
of the flat FRW model which would yield the same luminosity distances by inverting the FRW relation 



D L {z) = {l + z) 



dz' 



(2.10) 



to find 



H(z) 



d_ f D L {z) 
dz \ 1 + z 



(2-11) 



1 Note that the mass function M(r) which appears in Bondi Il8l . which we denote by Mg(r), is related to Celerier's M(r) by M' B (r) 
M'(r)/y/l + 2E[r), and so our radial coordinate specialization in Bondi's notation is M' B (r) = iirr 2 p/ \Jl — 2k(r)r 2 . 

2 More generally, an observer might fit data on Dl(z) to FRW models with arbitrary spatial curvature, including flat ones. 
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We can then calculate the associated deceleration parameter 



«(*) 



1 



H(z) 



dH{z) 
dz 



and the effective equation of state parameter 



w e g(z) 



?(*) 



2(1 + z) d 
3 dz 



In 



(1 + z) 3 / 2 



(2.12) 



(2.13) 



If we know to(r) and E(r), then we can find R{r, t) very simply by using the appropriate solution above, chosen from 
Eqs. I|2.2|) . (|2.3|) . and I2.4|l . We then solve the differential equations (|2.7I) and l|2.8|l to find t(z) and r(z), starting 
from the initial conditions r = and t — t a . We insert these t(z) and r(z) into the right hand side of Eq. (|2.9I) to 
obtain Dl(z), and then use Eqs. (|2.11l) and (|2.12l) to find H(z) and q(z). We will use this procedure later in this 
section with a class of models as an illustrative example. 



B. The Weak Singularity at r = 



There have been many claims that there exist LTB cosmological models in which go = q(z = 0) < [Tgt l20t l2lt 12^. 
l23.l2Sl23.l27| . For example, Iguchi et al. (2?J look at two different classes of LTB models: (i) models with k(r) = 
and a pure "BigBang time inhomogeneity" and (ii) models with to(r) = and a pure "curvature inhomogeneity" . 
In either case, they try to reproduce the luminosity distance function of a flat FRW universe with a matter density 
&m = 0.3 and a cosmological constant density f^A = 0.7, namely 



Dl(z) 



1 



dz' 



n M (i + z'y + n A 



(2.14) 



They appear to be successful up until they find R' < or R! < at a redshift z ~ 1 (we will discuss these pathologies 
in the next section). Thus, they appear to successfully find models where go < 0. 

On the other hand, the local expansions of Flanagan Q and of Hirata and Seljak show that qo is constrained to 
be positive for arbitrary inhomogeneous dust-dominated cosmologies that are not necessarily spherically-symmetric. 
In particular, Flanagan expands the luminosity distance as 



D L = A(6, 4>)z + B{9, </>)z 2 + 0{z 3 ) 



(2.15) 



where 9 and <f) are spherical polar coordinates as measured in the local Lorentz frame of the observer. He then defines 
the central deceleration parameter as 



qo 



where angle brackets denote averages over 9 and 
assuming that the pressure is zero, he finds 



1 - 2H 2 (A 
0, and Hq - 



- 3 B) 



(2.16) 

Using local Taylor series expansions and 



Qo 



47T 

3Hl 



P + 



3H* 



a/3 



a/3 



(2.17) 



where a a p and uj a p are the shear and vorticity tensors. The first term of this expression is obviously positive, and the 
terms in the brackets vanish in LTB models by spherical symmetry. Thus there is a contradiction: general theorems 
prove that go is positive in these inhomogeneous models, whereas the analysis of specific examples appears to show 
that it is possible to construct models in which go can be negative. Here we present the resolution of this contradiction, 
that there exists a weak local singularity which is excluded at the start from the computations of Flanagan and Hirata 
and Seljak, but which is present in models giving go < 0. We will show that the exclusion of this singularity inevitably 
leads to models with a positive go- 

We expand the density (|2.5|) to second order in r as 



p(r, t) = po(t) + P i{t)r + p 2 (t)r 2 + O (r 3 ) 



(2.18) 



The weak singularity occurs when pi(t) is nonzero, in which case the gravitational field is singular since OIZ — > 00 as 
r — > 0, where 1Z is the Ricci scalar. In other words, second derivatives of the density diverge at the origin, independent 
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of where observers may be located. This is true both in flat spacetime and in the curved LTB metric when we have 
a density profile of the form (|2.18l) . The singularity is weak according to the classification scheme of the literature 
on general relativity j2^| . This singularity is excluded from the start in the analyses of Flanagan [3] and Hirata and 
Seljak ^(j which assume that the metric is smooth. 

We now determine the conditions for a weak singularity to occur. We define the variable 

a(r,f) = ^M ; (2.19) 
r 

this is analogous to the FRW scale factor a(t), in the sense the metric takes the form 

ds 2 = -di 2 + a 2 (r, t) { [i + ^'MVaM)] 2 ^ + r 2 (d9 2 + sin 2 9d<f> 2 ) } . (2.20) 

1 — k{r)r z v ' I 

We expand this function as 

a(r,t) =ao(t) + cn(t)r + a 2 (t)r 2 + O (r 3 ) . (2.21) 
Comparing this to the formula (|2.2|l for R(r,t), we find for the zeroth order expansion coefficient 

a (t) = [6irGp (t o )} 1/3 t 2 / 3 . (2.22) 

We define Hq = cio(t )/ao(t ), and our choice of units above imply ao(t Q ) = 1. Using Eqs. (|2.19|l and l|2.21|l in the 
expression ()2.5(l for the density gives 



Pa (to) _ 4 /Qq (t ) oi (t) 
al(t) a 3 (t) 



p(r , t ) = ^1 - ^ n :°sZ r + O (r 2 ) . (2.23) 



Since ao(t) ^ by Eq. I|2.22|l . we see that having a non-singular model requires ai(t) = 0, or equivalently R"(r = 
0, t) = 0. 

It is straightforward to see that if a\ = 0, then qo > 0, and that if ai(t) ^ 0, then qq may be positive or negative. 
Note that the observer's measurement of qo from the data does not depend on the observer's prior assumptions about 
spatial curvature, and so the following analysis of qo is sufficiently general and applies for arbitrary k(r). If a±(t) = 0, 
then the angular size distance is R(r, t) — rao(t) + r 3 a 2 (t) + O (r 4 ) , where r and t are evaluated along the path followed 
by a radially directed light ray. Evaluating the redshift for such a ray gives to lowest order z = Hqt + 0(r 2 ). Thus, 
the angular size distance is unaffected by density gradients up to terms of order z 3 . In other words, the standard 
expansion of the angular size distance R = Da to order z 2 , 

H D A (z) = z - iz 2 (3 + q ) + 0(z 3 ) , (2.24) 

is completely determined by the evolution of the uniformly dense core region of the expanding spherically symmetric 
model, where the density is po(t) — po(t ) / a^t) from Eqs. I|2.6I) and 12.22fl . which is the density of dust expanding 
with scale factor ao(t). Therefore, the effective values of qo for such a model must lie in the same range as are found 
for exactly uniform, dust dominated FRW models: q > 0. 

We can gain further physical insight into the behavior of LTB models near r = by expanding the field equations 
in r, assuming (see Eq. (|2.21|l l 

a(r, t) = oo (t) + a n (t)r n + ... , (2.25) 

and correspondingly 

k(r) =k + k n r n + ... ; (2.26) 

we show in Appendix 1X1 that a±(t) = corresponds to having k\ = via a direct analysis of the LTB solutions. 
Thus, for non-singular models, n — 2 is the leading order correction to strict homogeneity near the center. The field 
equations are given in Bondi |18|. and in our notation his Eq. (24) is 

1 ( dR(r,t)\ 2 AnGpo (t a ) r 3 _ 1 



dt J 3R(r,t) 2 



= -zKry . (2.27) 
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Substituting R(r, t) = ra(r, t) we find 



da(r,t)\ 2 _8nG P0 (t o ) 



dt J 3a(r, t) 

Using the expansions i|2.25[) and (|2.26() and equating like powers of r, we find 

a _ 8irG Po (t ) _ H$n 2 

a Q — k = h ii (L — il ) 

da a 
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2d a„ = — -5 fc« = 2 fc ™ ■ (2.29) 



The first of Eqs. (|2.29|) is exactly the same as the Friedmann equation for the scale factor ao(t) in a universe with 
arbitrary spatial curvature, subject to the single physical requirement f^o > 0. To solve the second equation, notice 
that do = — iJQf2o/2aQ, so rewrite it as 

a a n - a a n = a — — = - — , (2.30) 
at \ eta J 2 



which has the solution 



L Jo a 0\ l ) 



where C is a constant. Let us define S n (t) = a n (t)/ao(t); then Eq. I|2.31|l becomes 



lo H 2 (t')a^(t') 2 7 ao(0) iP(a )ag 




where 

a_ /n. i n 

(2.33) 

Comparing with results in Peebles [2^|, we see that the first term of Eq. (|2.32|) is just the shrinking mode of linear 
theory, and the second is the growing mode. The amplitude C of the shrinking mode is related to the bang time 
function by to(r) — — Cr n + . . ., and the growing mode amplitude k n is related to the lowest order energy perturbation 
k n r n . Note that this approximate solution holds for a n r n << ao(t), i.e. for < r n << 1/S n . 

We have shown that for n > 2, the central value of qo is greater than or equal to zero. We now compare with the 
mildly singular case with n = 1. The evolutions of ao(t) and ai(t) are governed by Eqs. (|2.29|) and (|2.31l) . For this 
case, the low z expansion of Da{z) depends on ai(t ) and fci, and the effective value of qo near the origin becomes, 
from Eqs. ijO^I . tfFM . and QTZty . 

1 2ai(t ) . ai(*o) l n ai(*o) /„ , fa 



, = :S] -^™ + ^ = :fl -^ 2 + - (2.34) 
y ° 2 ° Ho 2 ° H \ 2 ) 2H§ 1 7 

This is no longer constrained to be positive. 

In Appendix 1X1 we show directly from the solutions to the Einstein equations, Eqs. I|2.2II . (|2.3|l . and l|2.4|) . that we 
must have both £q(0) = and fc'(0) = in order to have a non-singular model, and if this is true then q cannot be 
negative. In Appendix we use this to show that the models of Iguchi et al. have weak central singularities. 



C. Achieving a Negative Apparent Deceleration Parameter at Nonzero Redshifts 

Although models that have been previously analyzed contain central singularities, it is still possible to construct 
LTB models without such a singularity for which the effective deceleration parameter q(z), as defined in Eq. I|2.12|l . 
is negative for some nonzero redshifts. Here we explore a class of zero energy LTB models with a bang time function 
to(r) that is quadratic near r = 0, and therefore non-singular there. 
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In a zero energy LTB model, we have 

dt = -R'(r,t)dr (2.35) 

and therefore we can get the equation for t(r) along light rays that we observe from supernovae. Also, z is a function 
of r via Eq. Q2.8|) . specialized to k(r) = 0, and we get z as a function of r only by using our solution for t(r) along the 
rays. The bang time function is chosen such that it will (i) approach a constant for large r, so as to have a uniform 
density for large redshifts, and (ii) have no terms linear in r, so as to avoid a singularity at the center. Thus we 
integrate Eqs. (|2.8|) and (|2.35(1 with the bang time function choice 

to(r) = - 2 ^ cr 2 2 n2 (2.36) 

where A and D are dimensionless parameters, r c = [6nGpo(t )] 1 , and we choose units where r c = 1. We choose 
the initial conditions at the center, t(r — 0) = 1 and z(r = 0) = 0, and we integrate from the center outward. 

Figure n displays results for the effective q(z) that we calculate from the above model using Eqs. I|2.9|) . I|2.11|l . 
and (£E!) for various values of A and D, namely (A, D) = (0.094,0.14), (0.20,0.29), (0.46,0.62), (0.75,0.91), and 
(1.0, 1.2). We choose values of A and D for which the minimum value of q(z) that is attained is approximately — 1. As 
we can see, although all the models are forced to have q(z = 0) = 1/2, it is nevertheless possible for the deceleration 
parameter to become negative at nonzero redshifts, as we find a region of q(z) < for z < 1. 



4 







-1 



x=n ns4, n=n 14 j 




'. A=0.20, D = 0.89 j 




. _ A=0.46, D = 0.62 / 




. ._ X=0.75, D=0.91 / ' 




X=1.0, D=1.3 / 




/ ' 

// 
/ / 
1/ 

If' 

i 


s 

/ ■ 


mi ' // 





I V -I ~-T I I I I I I I I I I I I 

0.5 1 1.5 



FIG. 1: The effective deceleration parameter q versus redshift z for several quadratic bang time models (12.361 
which have a minimum q of approximately negative one. Plotted here are the data for models with (A, D) = 
(0.094,0.14), (0.20,0.29), (0.46,0.62), (0.75, 0.91), and (1.0, 1.2). 

In order to reproduce the current luminosity distance data, we want q(z) to quickly fall to from q(0) = 1/2 to 
q{z) f» — 1 and then stay at that value until a redshift z ~ 1. In Fig. [2 we plot several quantities that encapsulate 
some of the characteristics of the functions q(z), which are useful for assesing the feasibility of reproducing luminosity 
distance data. We define Az neg to be the width, in redshift, of the region where q is negative, and Az q< ~i to be the 
width of the region where q is below negative one. We also found that the large redshift behavior is unstable in these 
models: q blows up as we eventually approach the initial singularity. As an approximate measure of the location of 
this divergence, we define z max to be the redshift at which q(z) exceeds 3. Ideally, we want Az neg ~ 1, Az g< _i = 0, 
and z max — > oo. From Fig. |3 it does not appear as though this model can reproduce the data well, although it is 
conceivable that one could construct a model which gives more realistic results. We see that by increasing D, we also 
increase the size of the region with negative q(z), which makes the model more phenomenologically viable; however, 
by increasing D, we also decrease z max and thus make the model less physically reasonable. 
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FIG. 2: Several measures of the feasibility of quadratic bang time models, plotted versus D for A = 1, 0.751, 0.589, and 0.455. 
From top to bottom, we have plotted Az neg , Az 9 <_i, and z max , all versus D. 



III. THE INVERSE PROBLEM 



Unfortunately, it is highly unlikely that one could guess a bang time function fo(r) that would yield the experi- 
mentally measured luminosity distance Dl(z). A better approach would be to solve the inverse problem: given the 
appropriate D^(z), work backwards to try to find the corresponding tnfr), which may or may not exist. This approach 
has been taken before, but without avoiding the central singularity |2jj, as is shown in Appendix^! Models based on 
selected D^{z) generally break down at z *~ 1 upon encountering some pathology. We explore and clarify the possible 
pathologies below. 



A. General Properties 



In the LTB metric, the angular size distance is given by 

D A (r, t) = R(r, t) = rT 2 {r, t) , (3.1) 

where T = [t — toir)] 1 ^ 3 . Here we have specialized to units where 6TrGpo(t ) = 1. We also define the equivalent FRW 
radial coordinate to be 

r FRW (z) = (l + z)D A {z) , (3.2) 

in terms of which we have 

rT 2 (l + z)=r FR w(z) = ^^ . (3.3) 

1 + z 

Suppose we are given rpRw^)? and therefore Dl(z) and Da(z), and from this we wish to find the corresponding zero 
energy LTB model. 

The equations defining our flat LTB model with bang time function may be written in the form 



dT 1 dt ( 2r 1 
~d7 ~ ~3 + ~dr V9T3 ~ ^2 

and 

1 dz ( 2 2r dt 



(l + z)dr \ST 9T 4 dr 



(3.4) 



(3.5) 
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Multiply Eq. (EU> by 2/T and then add to Eq. JSTSJ) to find 

dt 3T d 



[T 2 (l 



dz 2(1 + z) (r/T- 1) dz 
we can also combine Eqs. (|3.4|) and (|3.5|1 such that we eliminate dtg/dr altogether to find 



r dT 



T dz 



^ + 7^-1 K" = 



dr 



1 



1 + 2 



Defining X = T 2 (l + z), these equations can recast into 



1 dX 
X~dz~ 



and 



dt 



(r-FRwy/TTI/X 3 / 2 - 1) 
r F Rw(3/2 - r FR wv / m/X3/ 2 ) 

3X 3 / 2 



3T 



3X 3 / 2 



drFRw 



2(1 



\3/2 



dz 2r FRW (l + z) 3 / 2 (3/2 - r F RwVTT7/X 3 / 2 ) 



3X 3 / 2 
2(1 + z) 3 / 2 



dz 



dr FRW 

dz 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



In the spatially flat, dust-dominated FRW model, X — 1 and rpRW^) — 3[1 — 1 / V 1 + z] . 

Given fpj^fz), Eq. I|3.8|) is a first order ordinary differential equation for X{z). It becomes singular when 



7-FRw(z)Vl + f 



for a flat, dust-dominated FRW model, this occurs when z = 5/4. Solutions z 
critical points of differential equation l|3.8|) . Near the critical point, 



(3.10) 

z cr ;t of Eq. (|3.1U|) . if these exist, are 



1 dX 
X~dz~ 



and 



dtp 

dz 



2(3/2- VTTz/x 3 / 2 ) 

?~FRw(z) 

(1 + z)(3/2 - r FRW \/Tn/X 3 / 2 ) 



1 + z 



d In r F Rw 
dz 



d In tfrw 



1 



dz 



(3.11) 



(3.12) 



Transcritical solutions, which are non-singular at the critical point, are possible provided that (l + z)dlnrpRw/dz = 1 
at the critical point. We discuss these solutions in more detail below. Clearly, the spatially-flat, dust-dominated 
FRW model is one special transcritical solution. For a general choice of rFR\v(-z)j however, the conditions for passing 
smoothly through the critical point will not be generically satisfied, and both d In X/dz and dtp/dz will diverge there. 
This suggests that a flat LTB model with a bang time function can only mimic a generic tfrw(z) up to some limiting 
redshift below z cr i t , where 



FRW 



(Zcrit) 



1 



d In r F Rw 



1 + Zcrit 



dz 



. 



(3.13) 



We shall argue below that only the special class of transcritical solutions can extend to infinite redshift. 
For exploring characteristics of the solutions, it proves useful to define the new variable 



V = 1 



2r F Rw\/l + z 
3X 3 / 2 



(3.14) 



Substituting Eq. 13.14|l into Eq. (|3.8(l gives, after some algebra, 



dV 

dz 



(1-4^ + V^ 2 ) (l-V0 2 dlnr F R W 



2V(l + z) 
(1 + V 2 ) _ 
2V 1 ■ 



2V 
d In tfrw 
dz 



dz 



1 + z 



din r FRW 
dz 



(3.15) 
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For a flat, dust-dominated FRW model, Vfrw = 3— 2V1 + z, and substituting this V[z) into the right hand side of Eq. 
(|3.15|) yields (LVfrw I dz = — l/\/i + z. Near z = 0, we have seen that flat LTB models resemble flat, dust- dominated 
FRW models, so 

dlnr FWW /dz = z^l - -(1 + q{0))z + ...} = z~ l (l - 3z/4 + ...), (3.16) 



2 



and therefore 



V(z)*l-z+--- + ... (3.17) 

for The first term in the small-z expansion of V(z) that can deviate from Eq. I|3.17|l is of order z A . 

At sufficiently small z, we expect V(z) to decrease. There are three possible classes of solutions to Eq. I|3.15|l : (i) 
solutions that decrease from V(0) = 1 to some constant Vx> < 1 as z — > oo, without crossing the critical point at 
V = 0; (ii) solutions that decrease until a redshift z = zq < z cr j t , where they terminate; and (iii) transcritical solutions 
that pass through the critical point smoothly. We examine these three classes in turn. In our considerations, we keep 
rvBw(z) general, with the provisos that the model tends to q = 1/2 at both z — > and z — > oo. The former is dictated 
by the character of LTB models free of central singularities, whereas the latter must be true of any phcnomcnologically 
viable model. In particular, then, we assume that H w §^frw(1 z) 3 ^ 2 at large z, where ^frw < 1- Therefore 
fFRw(z) — > f*FRW,oo as z — > oo, where rFRw,oo is a constant. 

Consider first solutions that decrease toward Voo asymptotically. At large values of z, Eq. (|3.15(l becomes 

dV 1-Woo + Vl (1- V oc ) 2 dlnr FRW 



dz 2Voo(l + z) 2Voo dz 

1 - 4Voo + Vl 3(1 - Voc 



2Voo(l + 2) 4V 00 ^ w r F Rw,oo(l + ^ 3 / 2 



(3.18) 



The first term on the right hand side of Eq. (|3.18|) is negative as long as V^ > V = 2 — \/3 « 0.27, and it dominates 
the second term. But V(z) ~ — ln(l + z) in that case, and this diverges. Thus, we can only have Voo = V). In that 
case, we let V(z) = Vq + u(z) at large z, and find 



dz 

which has the general solution 



du w ^3 3(V3-1) 2 (31Q) 



(2 - V3)(l + «) 4(2 - >/3)n^rpHw,oo(l + ^) 3/2 



«(z) = 3( ^- 1)2 , (3.20) 

where C is a constant. Although u(z) —* as z — > oo, it approaches zero from below, not from above, which 
contradicts our basic assumption. Thus, solutions that simply decrease toward constant V(z) > asymptotically do 
not exist. Conceivably, there can be solutions that decrease to a minimum and then increase toward Vq asymptotically. 
For these, however, V(z) will be double valued. It then follows that X = T 2 (l + z) must be double valued, since 
rpRwVl + z is monotonically increasing, and such behavior could be pathological. More generally, we shall see below 
that solutions that avoid V = must terminate at z = z C rit in order to avoid other physical pathologies. Thus, a 
solution "on track" to a minimum value V m i n > 0, and on to Vq asymptotically, might end at finite redshift. 

Next, consider solutions that reach V = at z = zq < z C rit and end there. Near V = 0, Eq. H3.15(l is approximately 



dV _ K FRW (z ) 

dz~ 2V ' [6 - 2L) 

where i^FRw( z ) = (1 + z)^ 1 — dhxr F -gsff /dz, as in Eq. I|3.13|l . Note that since zq < z cr ;t, -Kfrw^o) < as well. The 
solution to Eq. H3.21fl is V(z) ~ \/—Kfkw(zo)(zo — z), which terminates at z — zq. 
In order to reach the critical point, we must have 

K FRW (z) < (3.22) 
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all the way up to the critical point, with equality holding at z = z cr jt for the transcritical solution. For a transcritical 
solution to exist, we must be able to expand 



K Fnw (z) = QAz + 0(Az 2 



(3.23) 



near the critical redshift, z cr ; t , where Az = z — z CT it an d Q > 0. For a flat, dust-filled FRW model, we have z cr ;t = 5/4 
and Q = 8/27. Using this linear approximation, we find from Eq. (|3.15l) that V — kAz + 0(Az 2 ), where the slope 
k < is the solution to 



k 



k _ Q 

1 + z crit 2 







That is, we need the negative root 



k = 



1 



2(1 + Zcrit) 




(3.24) 



(3.25) 



For a flat, dust-filled FRW model, we find k = —2/3. This is clearly a transcritical solution. 

We can turn the above analysis into a test of whether a candidate for rp-nw( z ) that agrees with observations 
can be represented by a transcritical, zero energy LTB model. First, for the candidate model, it is possible to find 
z cr j t and Q algebraically; we can find z cr it using Eq. I|3.13|l . by requiring that i^FRw( z crit) = 0, and we can find 
Q = cLKfbw /dz\ z=ZcTit , cf. Eq. 13.23fl . Next, find k given Q and z cr ;t from Eq. 13.2511 and use this value of k to 
integrate Eq. Ij3.15|l back toward z = 0. If the solution satisfies Eq. (|3.17(l as z — > 0, then it is an acceptable 
transcritical solution. 

There are other disasters that may befall the solution for general f"FRw(z)j and some of these may even afflict 
transcritical solutions. Eq. (|3.9|) may be rewritten as 



dtp 

dz 



3(1 + *)(1-V) 



1 



1 



d In rpRw \ 



1 



y \i + z dz j (i + z)(i-v)_ 

As we have noted before, dtp/dz diverges at V = for generic j-frw(z)j but for transcritical solutions, 

dt 2r FRW (z cri t) / Q , 1 



dz 



3(1 + Z crit ) \k 1 + Z crit 



0{Az) 



(3.26) 



(3.27) 



near the critical point, which is finite, so this potential disaster is avoided. In particular, for a flat, dust-filled FRW 
model with Q = 8/27 and k = —2/3 at z cr it = 5/4, we see that Q/k + 1/(1 + z C rit) = 0> which is consistent with 
to(z) = for all redshifts. 

We must check for two other possible disasters, for solutions that are transcritical or not. As mentioned in the 
previous section, physical regions in any solution must have a positive, finite R' — dR/dr and dr/dz. We find 



dR 



( 2rpRw \ 



2/3 



(1 + z)d\ur-p-BSN / dz — 1 



dr 
dz 



2rpRw 



3(1- F)(l + z) 4 



1/3 



(1 - V)[(l + z)d]nr F Rw/dz -1] + 2V 
1 - V 



2V 



(1 + z) 1 

dz 



1 



(3.28) 



Note that we have R 1 oc (dr / dz)^ 1 , and thus a finite, positive dr/dz implies that R' > 0. These equations are 
evaluated along the path of a light ray directed radially inward. The requirement dR/dr > along the light ray is a 
necessary but not sufficient condition for an acceptable model. The more general requirement is that dR/dr > at all 
(r, t), a global condition that is much harder to satisfy; but in general, from Eq. I|2.2I) . this will be satisfied for models 
with to(r) decreasing monotonically. From the first of Eqs. (|3.28|l . we note that dR/dr — > at z = z cr it for solutions 
that are not transcritical. Solutions that terminate at zq < z cr i t would not encounter this pathology. Solutions of 
the first type described above, which decrease from V(0) = 1 but do not cross 1^ = 0, would end at z = z c rit- F° r a 
transcritical solution, 



dR _ 1/3 / 2r FRW (z crit ) \ 2/3 
dr \ 1 + 2 C iit / 



0(1 + ^crit) 
Q(l + Zcrit) - 2fc 



> 



(3.29) 
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at the critical point. Transcritical solutions therefore propagate right through the critical point with a positive, finite 
dR/dr. From the second of Eqs. (|3.28|) . we see that dr/dz diverges for solutions that terminate at V — and z = z . 
For transcritical solutions 

f-r + ^ — )>° ( 3 - 3 °) 

V k 1 + Zcrit/ 

at the critical point. Beyond the critical point, transcritical solutions have V < 0, and for reasonable rFRw(^) with 
decreasing (1 + z)dlnrpRw / 'dz , it seems likely that dr/dz remains positive. 

From these general considerations, we conclude that zero energy LTB models can only mimic a given, generic 
^frw(^) _ arranged, for example, to fit observations of Type la supernovae - for < z < zq < z cr ;t, where z cr it 
is the solution to Eq. i|3.13[l . There can be exceptional, transcritical models that extend to infinite z without any 
mathematical or physical pathologies. However, transcritical models are highly constrained mathematically, and may 
not exist for choices of rFRw(z) that conform to phenomenological requirements. The flat, dust-dominated FRW 
model is one transcritical solution, but it is ruled out by observations. 



dr 

dz 



^FRW^crit) 



3(1 + Z CT it) 



B. Manufacturing Transcritical Solutions 



To manufacture transcritical solutions, we will specify V(f), where f(z) = r-p^f(z)y/l + z, and find an equation for 
f{z). From Eq. (|3.15(l we find 

dr 3-10V + 3V 2 

dz 2{l + z)[2VV' + (l-V) 2 /f] [ ' 

where V'(f) = dV(f)/df. As long as W — > near V = 0, Eq. I|3.31|l satisfies the transcriticality condition -Kfrw = 
when V = 0. As an example, suppose we assume that V = 1 — kf. Then we find f = (2/k)(y/l + z — 1) and we 
must choose k = 2/3 in order to have the proper behavior at small z. This solution is simply equivalent to the flat, 
dust-filled FRW solution. 

Superficially, the prescription is simple: specify a V(f), make sure that VV — ► when V = 0, and then find the 
corresponding f(z) by integrating Eq. (|3.31|) . However, we know that acceptable solutions must have dr/dz > and 
finite; these conditions are not so easy to guarantee. 

Let us assume V(f) = 1 - §r/(f); then Eq. (|3~3H becomes 



3' 

dr (r/-l)(r/ + 3) 



dz 2{l + z)[f{ff - !)-{!- 2ff/2,)ff] 



(3.32) 



The numerator of Eq. (|3.32|l is zero when ff(f) = 1, or V = 1/3. If the denominator of Eq. (|3.31(l is nonzero at this 
point, then dr/dz goes to zero, and changes sign upon crossing it. Thus, we also want the denominator to vanish 
for an acceptable solution. In other words, V = 1/3 must be a critical point of Eq. (|3.31|l : we have only succeeded 
in hiding the critical nature of the problem, rather than eliminating it. Eq. (|3.32|l shows that to pass through this 
critical point we must require that ff = when ff = 1. Clearly, the spatially flat, dust-filled FRW model, for which 
/(f) = 1, is one possibility. 

It is also possible that the denominator of Eq. (|3.32l) vanishes, so dr/dz — » oo before ff — > 1. This happens when 

- wBwm ■ ( " 3) 

If /' < at small values of f, it is possible that infinite dr/dz occurs before ff — > 1. Since we also want /(0) = 1 for 
nonsingular models, and /(oo) = constant for models that approximate a flat FRW model with < 1 at sufficiently 
large redshift, we have several requirements on f(f) that must be met simultaneously for a model that is acceptable 
mathematically. Moreover, physically acceptable models must also have dR/dr > and R' > 0. Only a subset of 
such models - if any - will also be acceptable phenomenologically. 

To examine the phenomenological properties of a candidate transcritical solution, first define the effective Hubble 
parameter iJ c ff(z) via 

d ^L = 1 . (3.34) 
dz H e s{z) 
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It is straightforward to show that, normalizing so that H c s (0) 



h(z) 



H eS (z) 



(1 + 2) 



df 
dz 



(l + z)3/ 2 2 

We can also calculate the effective value of the equation of state paraemter 

2(1 + 2) din /i 

(f/ + 3)(f/-l) 



3[f(ff -!)-(!- 2ff/3)ff] 
3(f/ - 1) + (1 - 2ff/3)f 2 /' 



Wcff 



3[/(7 - 1) - (1 - f7)7'] 2 [3(f/ - 1) + (1 - |f/)f 2 /'] 



x /' 



4(f/) s 



+ 2(7) 2 - 8ff + 6 



2ff 



+ (rf) 2 



-f/"(f/-i)(7 + 3) U-^f 



2(7) 2 , 2f/ 
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3 



(3.35) 



(3.36) 



H c s(z) is the Hubble parameter that would be measured by an observer who assumes her observations are described 
by a spatially flat FRW model, and w e g(z) is the associated equation of state parameter. (A less dogmatic observer 
would allow for the possibility of spatial curvature.) Note that Eq. H3.35JI implies that for /' = 0, H e s/(l + z) 3 / 2 — 
f = constant. If we are interested in using LTB models to mimic a spatially flat FRW model with a mixture of dust 
plus cosmological constant, we shall want / — > v^T/u = \A ~~ at large redshift, where and Sly are the present 
density parameters in nonrclativistic matter and cosmological constant, respectively. Moreover, / — + 1 as z ~ * 0; thus 
/ must decrease from 1 to y/1 — f2y as redshift increases to mimic observations in a FRW model of this sort. 
For any choice of /(f) tailored to pass smoothly through both V = and f/(f) = 1, we can integrate Eq. 



to find a transcritical solution. However, such a solution still must pass the tests outlined in the previous section to 
extend to arbitrarily large redshifts. In terms of f and /(f), 



(1 



dtp 
dz 

dr 



3(h - /) 



3R' 



2(l + z) 1 /2 



2(1 + z) 5 / 2 f 2 hf(l 
/ 1 / 3 (3- 2hf) 
f + 2h- 2fhf 
3/V 3 (l - 2f//3) 
2 + f/h-2ff ' 



2f//3) 



(3.37) 



We reject any transcritical solution for which dto/dz is ever positive, or dR/dr or Rf change sign. 3 Moreover, even 
if a transcritical solution is found that possesses none of the pathologies discussed above, it may not conform to 
observational constraints. Thus, if there are any nonsingular, non-pathological LTB solutions that can also mimic the 
observations, they must be very exceptional indeed. 

Designing nonsingular, nonpathological transcritical solutions is a formidable challenge. Suppose that at small values 
of f, we expand /(f) = 1 + f n r n + . . .. Then we find that h w 1 + / n (l + n)f n + . . . and w c s — f n n(n + l)f™ _1 + . . . 



near the origin, where f 



Also 



(1 



dtp 
dz 
dR 
dr 



3R' 



2(1 



a/2 



3(A-/) „3 
2f ~ 2 
2{h-l) 



l + \{h-l) 



(3.38) 



3 For LTB with bang time perturbations only, (8R,/dr)t = [t — to(r)] _1 / 3 [t — to (f) — ^rdto/dr], which is only zero for a shell at coordinate 
radius r when t — to(r) = ^rdto(r)/dr. If dto(r)/dr < 0, this occurs before to(r) and is therefore irrelevant. As long as dto(z)/dz < 
and R' > along the light ray path, dto(r)/dr < and (8R/dr)t is never zero. 
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We wish to tailor /(f) to maintain positive values of both dR/dr and R', but already near the origin dR/dr and R' 
deviate from their fiat, dust-filled FRW relationships in opposite senses. Notice that to avoid the weak singularity 
near the origin, we need to have n > 2. Moreover, we want to make sure that to(z) is monotonically decreasing to 
avoid shell crossing. Near the origin, decreasing to(z) implies /„ < 0. 

To illustrate how difficult it is to manufacture non-pathological transcritical solutions from Eq. (|3.32|l . we have 
considered 



m = i 



1+K 



(3.39) 



The model has four parameters: A, fi, n, and p; the remaining parameter f 2 will be determined in terms of these 
four. To be sure that h and w e g are finite near ff = 1, we want either p = 2 or p > 3. Expanding near the origin, we 
find 



l + K(h/h) n _ Kp /V "" 
h V ^2 



1 



1 



requiring /(0) = 1 inplies 



K 



f x - 1 



and so we can rewrite the expansion in the form /(f) = 1 + /«T n 



with 



h \f 2 J If™ 



p{h - 1) 
' i 



ri 



A" 



(3.40) 



(3.41) 



(3.42) 



Thus, we want f i — 1 > and therefore K > for /„ < and real f 2 . At large values of f, /(f) — > ff x [l + (— 1) P A] . 
Thus, we expect models that can mimic decelerating FRW models successfully to have f f 1 [l + (— 1) P A"] < 1, suggesting 
either large f i or odd p, or both. Empirically, we have been unable to find any non-pathological models based on Eq. 
()3.39(l with these properties. 

Figure |21 shows an example of a transcritical solution with K = 1, fi = 1.05, n = 3, and p — 2. Although the 
figure only displays z < 1000, we have integrated this model out to z = 10 6 to verify that it asymptotes smoothly to 
a high redshift FRW model, with constant t . The left panel shows »TRw(^)/ r FRw( z ) (dotted line), where ^^(z) 
is computed for a flat ACDM FRW reference cosmology with Qm = 0.27 and fl\ = 0.73, h(z) (short dashed line), 
and w c ff(z) (solid line). Since h — > 2/1.05 = 1.905 at high z for this model, inevitably there must be regions with 
w c ff > 0; this is in the range 0.75 < z < 3.9, with a peak value w c r ~ 2.13. There are two regions of negative w e s- 
(i) one at < z < 0.75, with minimum value u> c ff ~ —0.98; and (ii) an extensive regin at z > 3.9, with a minimum 
value w c s ~ —0.292. The right panel verifies the transcritical nature of the solution: it shows V (short dashed line), 
1 — (ilnrFRw/^ m (l + z) (dotted line), and to (solid line), and the long dashed line is at 0. The first two cross zero 
simultaneously, as they must for a transcritical solution, and at large rcdshifts, V is approximately proportional to 
\/l + z while 1 — dlnrFKw/rfhi(l + z) tends toward one. For this model, to(z) < at all z, and we have also verified 
that it decreases monotonically. In addition, we can verify that the model behaves as predicted at small redshifts: 
t (z) fa -0.342 3 and w cS m -2.7z 2 . 



Figure 0] compares the relative distance moduli 



: 5.01og 10 [r FRW (z)/4°j{ w (z)] 



(3.43) 



for models with K = 1 and (n,p) — (3,2) (solid line), (n,p) — (3,4) (dotted line), (n,p) — (2,2) (dashed line), and 
(n,p) = (2,4) (dash-dot line), with ?\ — 1.05 in the lower set of curves and f\ = 1.5 in the upper set; there is no 
solid line in the upper set for (n,p) = (3,2) because the model is pathological. For the lower set, luminous objects 
would appear systematically brighter than they would in the standard ACDM model. As f% is increased, a period of 
substantial acceleration is seen in the models below z ~ 1, leading to the systematic brightening relative to ACDM 
as seen in the upper set of curves. In either case, the luminosity differences would be easy to discern observationally. 

These few models illustrate several important qualitative points. First, it is possible to construct non-pathological 
transcritical solutions that also avoid any central weak singularities. Second, the model has a complicated "effective 
equation of state", including regions where w e g < 0, but also regions where w e g > 0. In this case, we found a range 
of values — 1 < w c g < 2. Finally, although it is possible to construct models that are well-behaved mathematically, 
these models do not generally conform to observational constraints. We have not, however, excluded the possibility 
that transcritical models in agreement with observations may exist. 
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FIG. 3: Results for a transcritical model based on Eq. I|3.39jl with K = 1, p = 2, n = 3, and fi = 1.05, out to z = 1000. 
The left panel shows tf-rw(z) /^^(z) (dotted), h(z) = H(z)/[Ho(l + z) 3 ^ 2 ] (short dashed), and w c g(z) (solid). The reference 
model corresponding to 4rw( z ) is t ne spatially flat ACDM model with SI m = 0.27. The right panel shows V (short dashed), 
1 — dlnrFRw/rfln(l + z) (dotted), and to (solid). 




FIG. 4: Distance moduli relative to spatially flat ACDM with Qm = 0.27 for models with K = 1 and (n,p) = (3,2) (solid), 
(n,p) — (3,4) (dotted), (n,p) = (2,2) (dashed), and (n,p) = (2,4) (dash-dot), with fi = 1.05 for the lower set of curves, and 
fx — 1.5 for the upper set. 



IV. CONCLUSIONS 



Observations of Type la supernovae imply that we live in an accelerating universe, if interpreted within the frame- 
work of a homogeneous and isotropic cosmological model 0, 0, ■ Some have tried to use the spherically symmetric 
LTB cosmological models to explain this seemingly anomalous data, as introducing a large degree of inhomogeneity 
can significantly distort the dependence of luminosity distance on redshift. We have shown that one must take care 
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in using these models, as they will contain a weak singularity at the symmetry center unless certain very restrictive 
conditions are met. Realistic LTB solutions require that the first derivative of the bang time function vanish at the 
center, t' (0) = 0, and also that fc'(0) = 0, where 2E(r) = —k(r)r 2 . Otherwise there are physical parameters, such as 
the density and Ricci scalar, which are not differentiable at the origin. 

We have also shown that any LTB models without a central singularity will necessarily have a positive central 
deceleration parameter go, and thus all previously considered LTB models with qo < are singular at the origin. 
However, it is still possible to obtain a negative effective deceleration parameter for nonzero redshifts, which we have 
shown using as an example the model with zero energy and with the bang time function i|2.36[) . that is quadratic 
at small r. These models have regions of apparent acceleration, where q(z) is negative. If our goal is to reproduce 
luminosity distance data with a non-singular LTB model, we can try to smooth out the center appropriately and 
tailor the model to fit the data at modest redshifts, say z > 0.01. This is not an easy task because there are other 
singular behaviors that generally occur when trying to represent a given luminosity distance function Dl(z) with a 
zero energy LTB model. 

Our detailed examination of the "inverse problem" elucidates how difficult it is to match zero energy LTB models 
to observed luminosity distance data. We have shown that the underlying differential equations generically become 
singular at a critical point. We have also shown that some exceptional choices of tfrw — Dl(z)/(1 + z) admit 
transcritical solutions which are smooth at the critical point z — z cr j t , and may extend to arbitrarily high redshift, 
given that they do not encounter other pathologies along the way. All other solutions terminate at some redshift 
z < z cr it- We have shown how transcritical solutions can be constructed via a simple procedure. Although these 
solutions show both enhanced deceleration, seen as regions with w e s(z) > 0, and acceleration, seen as regions with 
w e s(z) < 0, none that we have constructed explicitly conform to observations. Here we have only studied the effects of 
a bang time function, and did not consider the case of a non-zero energy function E(r). We expect generic solutions 
with E(r) 7^ to share the basic characteristics of the models studied here, namely the critical points and other 
singularities that we have discussed. However, we cannot say for sure that there are no transcritical and nonsingular 
solutions with non-zero to(r) and E(r) that agree with observational data on rFR.w (z). although it does not appear 
to be likely, as is evident from previous unsuccessful attempts to find such solutions |23 ■ 

Even if it were possible to reproduce determinations of Dl(z) from supernova data in a LTB model without 
dark energy, we would still be left with the task of matching all of the other cosmological data with such a model. 
First, the Wilkinson Microwave Anisotropy Probe is one of our most important sources of information about the 
Universe, via CMBR data; Alnes et al. |30j try to reproduce the first peak of the angular power spectrum with LTB 
models, and Schneider and Celerier [3l| claim to be able to account for the apparent anisotropy in the dipole and 
quadrupole moments with an off center observer. There are further constraints on inhomogeneous models from the 
kinetic Sunyaev-Zel'dovich effect, which constrains radial velocities relative to the CMB 32] . However, observations 
of large scale structure formation may be the most difficult to reconcile. These data strongly disfavor a currently dust 
dominated universe, as density perturbations would have grown too much without dark energy present to speed up 
the cosmic expansion rate and consequently retard the growth of fluctuations. 



In this appendix we show directly from the solutions Q2.2|) . (|2.3|) . and (|2.4|l of the Einstein equations that LTB 
models without central singularities must have positive go- The zero energy solution k(r) — has 



Thus, we see that the zero energy solution requires t' (0) = in order to have no central singularity. More generally, 
for the k(r) > solutions, we find 
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APPENDIX A: PROOF OF g > DIRECTLY FROM LTB SOLUTIONS 



ai (t) cx R"(r = 0, t) = - | [6nG Po (*o)] 1/3 *o(«)^ 1/3 ■ 



(Al) 




3F'(x ) 4 8ttGpo (t ) F fop) 

rz — £ ^ ■, n 



t' (0) 2^/h)F'(x a ) 



(A2) 



\fk~o 3fcp 
where k = k(r — 0) and we have defined the function F(x) by 

1 — cos u = F(u — sin u) . 



(A3) 
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Here Xo(t) — uq — sinito is the value of x at the center r = at time t: 



x (t) 



3k, 



3/2 



A-KGp (t Q ) 

Similarly, for the k(r) < solutions we find 

3G'(x ) 8nGpQ (t a ) G (x Q ) 



R"(0,t) = -fc'(O) 
where we define the function G(x) by 



V-ko 



-t - 



3fcQ 



-t'o(O) 2^k a G'(x ) 



coshu — 1 = G(sinhu — u) . 



(A4) 



(A5) 



(A6) 



Since the bracketed expressions in Eqs. (|A2|I and (|A5|) are functions of time, R"(r — Q,t) will vanish at arbitrary t 
only if t' (0) = and fc'(O) = 0, and only then can one avoid having a singularity at the symmetry center. 

These conditions, t' (0) — and k'(0) = 0, lead to the restriction that qo must be positive. Celerier j2|J expands 
the luminosity distance for small redshift and finds the second order coefficient to be 



D 



(2) 



1 


\d 2 D L - 


1 


R' 


2 


. dz2 . 


r=0 ^ 


R> 



R" R" 



R'R' 



R' 2 R'R' R' 2 



(A7) 



r=0 



where overdots again denote partial derivatives with respect to time. The deceleration parameter at r = is therefore 

R'R' R" R" ' 



qo = l- 2H D { ? ] = 



R' 2 R'R' R' 2 



(A8) 



If R"(Q,t) = to avoid a singularity, we find that the last two terms in the above expression are also zero, and we 
obtain 



Qo = 



R'R' 
R' 2 



J r=0 



Since R'(r,t) > to prevent shell crossing, and R 2 is obviously positive, we would need to have 

R'(0,t) = & {t) > 
in order to have a negative qo- For the k(r) = solution, 

\2irGpo (to) 1 1/3 



R'(0,t) 



t-^ < 



moreover, the k(r) > solution has 



and the k(r) < solution has 



R'(0,t) = 



R'(0,t) 



3ko V dt 



AnGpo (t ) ( du 



iirGpo (t Q ) ( duo" () 



< 



3fco \ dt 



(A9) 



(A10) 



(All) 



(A12) 



(A13) 



Therefore we can conclude that, in the absence of weak central singularities, all LTB solutions have positive qo since 
R'(0,t) is always negative. 
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APPENDIX B: MODELS OF IGUCHI, NAKAMURA AND NAKAO 



In this appendix we verify explicitly that the models with qo < studied by Iguchi et al. |27| contain weak 
singularities. For the first case in Iguchi et al., the pure bang time inhomogeneity, there will be no singularity if 
t' (0) = 0, as shown from Eq. (|A1|) . If we expand D L {z) for this FRW model in a power series around z = 0, we can 
compare this term by term to the expansion of the luminosity distance for a zero energy LTB model to find [2(j 



si 



A I 



1 + 5 



and 



0\ = -- 



t' (0) 29ff(Q) 5tff(0) 
a/3 2 4 a 2 /3 4 2 a 2 /3 

l*o(0) , 29^ 2 (0) , 5*ff(0) 



(Bl) 



+ 



2 a/3 2 8 a 2 /3 4 4 a 2 /3 ' 
where a = [67rG / 9o(t )] 1 / 3 and (3 = t 1 /^. Using the fact that Sl M + SIa = 1, we combine Eqs. (|B1|) and (|B2|) to find 

1 



(B2) 



rap SI A ■ 



(B3) 



A nonzero SIa requires that t' (0) is also nonzero, and hence there will be a singularity in such models. 

Iguchi et al. also look at models with to(r) = and positive E(r). By combining and rearranging Eqs. (6) and (39) 
from 261, we find that 



3ft a - 1 R'R! 



R'' 



R' 2 R'R' 



RT_ 
R' 2 



(B4) 



Plugging into this the negative k solution for R(r, t) and then setting r = 0, we can find an equation for fc'(0). Iguchi 
et al. make some simplifying definitions, wherein they set Hq = G = 1 and then write everything else as a function 
of a parameter Oq, which they vary between 0.1 and 1. They set fco = SIq — 1, po(t ) — 3f2o/87r, 



Uq = 111 




and 



t(r = 0) 



SIq (sinhwo — uq) 



(B5) 



(B6) 



2 (1-Sl f 2 ' 

where t(r) is evaluated along radially inward-moving light rays. Plugging these in and then solving for fc'(0) yields 



fc'(0) 



(1 - Sl f 2 



G 



(3J7a — 1) sinh 2 Uq (coshu — 1) + 2 (coshu — 1) 
3sinhwo — uq (coshwo + 2) 



(B7) 



where it is assumed that Sljs. = 0.7. This shows that fc'(0) is only zero if Slo — 1; we can see from their Fig. (4) 
that this corresponds to the uninteresting FRW dust solution fc(r) = for all r. All of the other choices for Slo will 
correspond to models with fc'(0) ^ and a central singularity. Therefore, all of the non-trivial models computed in 
Ref. [23| have weak singularities at r = 0. 
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